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INVESTIGATION  OF  UNDERGROUND  EXPLOSIONS  WITH  MODEL  TESTS 
Measurements  on  the  platform 


SUMMARY 

Model  tests  in  three  scales  were  undertaken  to  examine  if 
simple  scaling  lavs  can  be  used  to  investigate  underground 
explosions*  Charges  were  detonated  in  steel  tubes,  and  pressure 
and  arrival  times  were  measured  in  the  tube  and  cn  a platform. 
The  objective  of  the  tests  is  described  in  more  detail  in  (l), 
which  also  describes  the  results  from  the  measurements  in  the 
tube . 

This  report  describes  the  results  from  the  measurements  cn  the 
platform.  As  for  the  measurements  in  the  tube,  a deviation  from 
the  simple  scaling  laws  is  observed*  The  dependence  of  the 
pressure  on  tube  diameter,  charge  weight  and  vail  roughness  is 
examined* 

A detailed  description  of  the  experimental  set-up  is  given  in 
(2)  and  of  the  data  processing  in  (3)*  Reference  (3)  also  con- 
tains a list  of  the  measured  front  pressures,  maximum  pressures 
and  arrival  times* 


1 INTRODUCTION 


Model  tests  in  three  scales  were  undertaken  to  examine  if  simple  scaling 
lavs  are  valid  for  describing  the  blast  wave  originating  from  an  acciden- 
tal explosion  in  an  underground  ammunition  storage  site.  The  objective  cf 
the  tests  is  described  in  more  detail  in  (l).  The  models  consist  of  a 
tube,  simulating  a tunnel*  and  a platform,  simulating  flat  terrain  (l,2)* 
The  pressure-time  history  of  the  blast  and  the  arrival  time  of  its  shock 
front  were  measured  and  recorded  (2,3).  The  results  from  the  measurements 
in  the  tube  section  were  presented  and  discussed  in  (l)  In  this  report 
the  results  from  the  measurements  cn  the  platform  section  are  discussed* 
The  necessary  definitions  for  reading  this  report  are  also  found  in  (1). 

In  the  appendix  to  (3) * all  the  recorded  front  pressures , maximum  pres- 
sures and  arrival  times,  both  from  the  tube  section  and  from  the  platform 
section,  are  listed* 
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DIFFICULTIES  WITH  PRESSURE  TRANSDUCERS 


The  pressures  measured  on  the  platform  range  from  a few  tar  close  to  the 
tube  outlet  to  about  h'O  mbar  at  larger  distances  from  the  outlet.  It  is 
difficult  to  find  good  pressure  transducers  which  cover  this  range  of 
values , 


We  have  used  LC65  transducers  from  Celesco . The  uncertainties  in  the 
measurements  with  these  should  be  about  20%.  In  addition  to  this  uncer- 
tainty, there  are  oscillations  in  the  platform  which  may  disturb  the 
measurements. 


When  the  experiments  were  finished,  it  was  realised  that  the  calibration 
values  of  the  transducers  as  given  by  the  producer  were  unsatisfactory* 
Control  calibrations  performed  at  the  Norwegian  Defence  Research  Estab- 
lishment  showed  a deviation  of  about  100$  from  those  given  by  the  produ- 
cer. No  information  about  linearity  is  given  by  the  producer. 

As  a result  of  this  uncertainty  with  the  calibration,  the  discussion  in 
this  report  is  based  mainly  on  the  arrival  time  measurements. 

However,  for  some  of  the  pressure  transducers,  it  was  possible  to  estab- 
lish a confident  calibration  procedure.  As  a result  it  was  possible  to 
recalculate  the  front  pressure  at  the  measurement  stations  closest  to 
the  tube  outlet  in  CFG . 1 , end  for  the  case  of  zero  wall  roughness  in 
CFG.l  even  at  most  of  the  measurement  stations.  The  front  pressures  ob- 
tained in  this  way  are  used  with  some  confidence. 

In  the  listing  of  the  pressures  in  the  appendix  of  (3)  the  uncorrected 
pressure  values  are  given. 
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3 TUBE  OUTLET 


3 « -1  Scale-dependence  of  the  pressure  drop 

In  (l)  we  observed  that  the  front  pressures  in  the  tube  were  larger  in 
the  larger  scales  , and  that  the  difference  increased  with  distance  fi*om 
the  charge.  When  we  look  at  the  front  pressures  at  the  measurement  sta- 
tion 10  tube  diameters  from  the  tube  outlet  on  the  platform  (Table  3.1) , 
it  is  therefore  surprising  to  find  that  the  front  pressures  now  are 
smaller  in  the  larger  scales.  The  deviation  from  scaling  with  respect  to 
dimension  seems  to  change  direction  when  passing  through  the  tube  outlet. 
F_gure  3.1  shows  some  examples  of  how  the  front  pressure  falls  as  it  passes 
through  the  outlet  in  the  three  scales.  We  shall  study  in  greater  detail 
the  pressure  drop  of  the  shock  front  when  it  passes  through  the  tube  out- 
let. This  pressure  drop  is  due  to  the  fact  that  the  propagation  of  the 
shock  wave  changes  from  one- dimensional  to  three-dimensional. 

Since  on  the  platform  we  have  confidence  only  in  the  pressure  measure- 
ments in  CFG.l,  we  shall  restrict  ourselves  to  this  configuration. 


3 . 2 Examination  of  the  pressure  drop 

We  have  calculated 


where  P'  is  the  front  pressure  at  position  n~52  in  the  tube  just  at  the 
outlet,  and  P' ' is  the  pressure  at  the  position  N-IO  tube  diameters  from 
the  outlet  on  the  platform.  P'  increases  with  increasing  scale,  while  P' ' 
decreases  with  increasing  scale. 

The  values  of  y are  listed  in  Table  3.2. 

To  study  the  deviation  of  y with  increasing  D,  the  values  of  y are  used 
to  estimate  a linear  equation 


Y 


a + bD 


(3.2) 


The  valuer  of  a.  and  b and  the  correlation  coefficient  R are  listed  in 
Table  3.3.  b is  positive,  indicating  that  the  pressure  drop  increases 
when  D increases.  It  is  also  seen  as  a weak  tendency  that  b decreases 
when  the  charge  group  q increases.  The  scale  dependence  of  the  pressure 
drop  is  thus  larger  for  smaller  charges  than  for  larger  charges. 

From  Table  3.2  we  also  see  that  y decreases  when  q increases,  and  that 
it  increases  with  increasing  large-scale  roughness. 


3 * 3 Conclusion 

The  large-scale  roughness  of  the  tube  walls  influences  the  pressure  drop 
in  a region  where  it  is  not  physically  present.  The  vail  roughness  must 
thus  influence  the  pressure  drop  through  its  earlier  impact  on  the  pres- 
sure-time history  of  the  blast  following  the  front  shock,  as  it  propagate 
through  the  tube.  It  must  be  concluded  from  this,  and  from  the  discussion 
in  section  3.2,  that  the  drop  in  front  pressure  through  the  tube  outlet 
depends  critically  on  the  pressure-time  history  following  the  shock,  and 
that  the  scale  dependence  of  the  pressure  drop  must  result  from  a scale 
dependence  of  this  pressure-time  history.  This  shows  that  if  one  wants  to 
get  a clear  picture  of  the  front  pressure  dependence  on  scale,  one  must 
examine  how  other  parameters  describing  the  blast,  such  as  impulse  and 
positive  duration,  depend  on  scale, 

We  have  not  analysed  our  pres sure /time  plots  sufficiently  to  give  any  con 
elusions  about  the  scale  dependence  of  impulse  or  positive  duration,  and 
shall  not  comment  on  which  physical  effects  are  determining  the  pressure 
drop.  As  will  be  seen  from  the  discussion  of  the  shock  propagation  on  the 
platform,  however,  the  scale  dependence  of  the  pressure  drop  in  the  tube 
outlet  does  not  seem  to  result  from  a scale  dependent  loss  of  energy  in 
this  region,  but  from  a scale  dependent  redistribution  of  energy  in  the 
blast . 
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CURVE  -FITTING  ON  PRESSURES  OBTAINED  FROM  TIME  MFASUREMEHTS 

Curve  fitting  procedure 


Based  on  measurements  of  arrival  times  a mean  velocity  of  the  shock  he- 
tween  two  succeeding  measurement  stations  are  calculated#  Yfe  let  this 
velocity  correspond  to  a position  midway  between  the  two  measurement 
stations.  The  front -pres sure  in  this  position  is  then  calculated  by  the 
Rankine-Hugoniot  equation.  The  curve  from 


P 


ov20' 


(4,1) 


where  N is  the  distance  from  the  tube  outlet,  measured  in  number  of  tube 
diameters,  is  fitted  to  the  values  of  the  front  pressures.  The  fitting 
parameters  P and  B represent  the  pressure  in  position  N ~ 20  and  the 
attenuation  of  the  front  pressure,  respectively*  Only  pressures  in  posi- 
tions H > 20  are  included  in  rhe  curve  fitting. 


The  values  of  P , B and  the  correlation  coefficient  R are  shown  in 
o 

Table  U.la  for  CFG.l  and  Table  l».lb  for  CFG. 2. 


Discussion  of 

It  is  seen  that  P^  increases  with  increasing  q,  but  that  the  increase 
is  smaller  than  a linear  increase.  The  curve  form 


P 


o 


bo 

a e 


(h.2) 


is  with  good  correlation  fitted  to  the  values  of  P . The  results  are 

o 

shown  in  Table  U.2a  for  CFG.l  and  Table  4,2b  for  CFG, 2. 


For  each  configuration  and  each  degree  of  the  large  scale  wall  roughness 
r^,  b seems  to  be  independent  of  D.  It  decreases,  however,  with  increas- 
ing rT • That  is,  while  the  increase  in  ? with  increasing  charge  group 
seems  to  bo  independent  of  scale,  it  is  smaller  for  larger  values  of  wall 
roughness . 


From  Table  h.2  it  is  seen  that  a,  on  the  average,  decreases  with  increas- 
ing rT . It  is  also  seen  that  a decreases  with  increasing  D,  and  this  do- 

L 
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crease  seems  to  be  larger  for  larger  values  of  rT . It  follows  that  P 

aJ  O 

decreases  with  increasing  r_  and  D,  and  .that  the  decrease  in  P is  larger 

JLi  O 

for  larger  values  of  This  is*  consistent  with  the  observations  in 
chapter  3* 

k . 3 Discussion  of  B 

From  Table  U.i  it  is  seen  that  E is  nearly  independent  of  q,  indicating 
that  the  attenuation  of  the  front  pressure  on  the  platform  is  nearly 
independent  of  a* 

The  equation 


B “ b]  + b2D  + b3rL  (h.h) 

is  by  linear  regression  fitted  to  the  values  of  B.  The  resulting  values 
of  b^}  and  bQ  and  the  corresponding  partial  correlation  coefficients 
are  given  in  Table  4,3a  for  CFG.l  and  Table  4.3b  for  CFG, 2.  All  the 
correlation  coefficients  are  significant  on  a level  e = 10 %, 

With  the  values  of  the  b-paraneters  inserted  in  Equation  (4.4),  we 
obtain 


B « 1.75  " 2.70  D - 5.39  rL 
for  CFG.l,  and 

B » 1.52  “ 2.01  D - 4.20  rL 


(4.5a) 


(4.5b) 


for  CFG. 2. 


We  tried  to  include  a term  b\ rrP  in  Equation  (4.4)  to  examine  the  inter- 

4 Lt 

action  between  roughness  and  scale.  The  resulting  values  of  b*  were  nega- 
live  both  for  CFG.l  and  CFG. 2,  but  the  correlation  coefficients  were  not 
significant . 


10 


h.3*l  Dependence  of  B on  scale 

The  first  thing  we  notice  is  that  B decreases  when  D increases;  that  is, 

the  attenuation  of  the  front  pressure  is  smaller  in  the  larger  scales. 

We  have  observed  that  P # the  front  pressure  in  position  N B 20,  decreases 

o 

with  increasing  scale.  The  situation  can  then  be  summarized  as  follows: 

In  the  tube,  the  front  pressure  is  larger  in  the  larger  scales.  When  the 
shock  wave  propagation  changes  from  one-dimensional  to  three-dimensional 
in  the  tube  outlet,  there  is  a large  pressure  drop,  and  this  pressure 
drop  is  so  much  larger  in  the  larger  scales  that  in  position  IT  = 20  on 
the  platform#  the  front  pressure  is  largest  in  the  smallest  scales.  As 
the  shock  propagates  on  the  platform#  however,  the  front  pressure  decrease 
more  in  the  smaller  scales,  until  it  eventually  becomes  smallest  in  the 
smaller  scales.  This  shows  that  the  scale  dependence  of  the  pressure  drop 
as  the  shock  front  passes  through  the  tube  outlet#  is  not  due  to  a scale 
dependent  energy  loss,  but  a scale  dependent  redistribution  of  energy* 

H.3«2  Dependence  of  B on  roughness 

We  see  that  B decreases  when  rT  increases*  The  effect  of  the  wall  rough - 
ness  in  the  tube  on  the  attenuation  of  the  front  pressure  on  the  plat- 
form must  be  through  its  effect  on  the  shape  of  the  pressure-time  profile 
of  the  blast  as  this  propagates  through  the  tube,  A small  attenuation 
corresponds  to  a flat  pressure-time  profile.  The  effect  of  the  roughness 
in  the  tube  must  be  to  flatten  the  pressure  time  profile  which  again  re- 
sults in  a smaller  attenuation  of  the  front  pressure  on  the  platform. 


k . k Results  from  front  pressure  measurement s 

As  mentioned  in  chapter  2 the  uncertainty  in  the  calibration  of  the 
pressure  gauges  is  generally  too  great  for  the  pressure  measurements  to 
be-  used  with  confidence.  For  r.  - 0 in  CFG.l,  however , the  necessary 
corrections  have  been  made,  and  the  resulting  values  of  the  front  pres- 
sures are  listed  in  Tables  U.Ha-c. 

Equation  ( t . 1)  is  fitted  to  these  front  pressure  values,  and  the  result 

ing  values  of  P and  B,  and  the  correlation  coefficients,  R,  are  also 
o 

listed  in  the  Tables  h.k. 


Equation  (<4.2)  is  fitted  to  the  obtained  values  for  P__ , and  the  results 
are 


P « 0.187  e 
o 


0.0030  q 


with  correlation  R “ 1.00 


for  scale  1:100, 


P *»  0.255  e 

o 


0.0022  q 


with  R ~ 0.99 


for  scale  2:100,  and 


P = 0.217  e 
o 


0.0025  q 


with  R ® 0.98 


for  scale  3:100. 


When  this  is  compared  with  the  corresponding  results  from  the  front 
pressures  obtained  from  arrival  time  measurements,  section  1.2,  Table  4.2 
it  is  seen  to  be  a fair  agreement  as  far  as  attenuation  is  concerned. 
However,  the  absolute  values  are  appreciably  lower. 


As  is  the  case  for  curve  fitting  based  on  front  pressures  obtained  from 
time  measurements,  B in  Equation  (4.1)  is  seen  to  be  nearly  independent 
of  q also  in  the  curve  fitting  based  on  the  measured  front  pressures. 

The  mean  value  of  B in  each  scale  in  the  latter  case  corresponds  reason- 
ably well  with  the  values  in  the  former  case. 


FITTING  OF  EMPIRICAL  EQUATIONS  TO  ARRIVAL  TIME  DATA 

Ir.  chapter  4 the  measured  arrival  times  were  used  to  calculate  the  front 
pressures,  and  a curve  fitting  procedure  was  performed  on  the  front  pres 
sure s . In  this  chapter  the  curve  fitting  will  be  performed  on  the  mea- 
sured arrival  times  directly. 


Comparison  between  c urve  forms 

Three  expressions  have  been  fitted  to  the  arrival  times  measured  in  posi 
tions  with  N > 20  on  the  platform.  The  expressions  are  the  following 
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B|  ( N- 20 ) + B2(N-20)2  + B3(lI-20)3 


(5.1) 


t-t 

D(H-20)  * aw(t-tQ)  + Aq  In  (l  + — “)  + Dq  (5-2) 

o 

t-t.  i 

D(N-20)  = a (t-t  ) + A {Infl  + — (5.3) 
0 0 o v B • 


t is  the  arrival  time  at  position  II  = 20.  a is  the 
o * o 

the  undisturbed  air  ahead  of  the  shock  front, 
are  fitting  parameters* 


speed  of  sound  in 

A > B and  D 
3 os  o o 


There  is  no  physical  reason  for  equation  (5*i);  it  is  only  a fitting 
equation#  Equations  (52)  and  (5*3)  have  at  least  the  property  that  the 
shock  velocity  approaches  the  velocity  of  sound,  as  time  approaches  in- 
finity# Expressions  similar  to  equations  (5*2)  and  (5<3)  vere , among 
others s,  discussed  for  spherical  explosions  in  free  air  in  a paper  by 
T K Groves  at  Suffield  Experimental  Station  (4). 


To  decide  which  of  these  three  expressions  that  give  the  best  fit  to  the 
measured  arrival  times,  a method  described  by  Dylewski  in  (5)  is  used. 


For  each  shot  the  curve  forms  are  fitted  to  the  observed  arrival  times, 
and  the  square  deviation  S is  calculated  for  all  the  three  curves.  Then 
the  measured  arrival  times  in  the  measurement  stations  nearest  to  the 
tube  cutlet  are  fitted  to  the  data  and  the  square  deviation  S*  is  calcu- 
lated for  the  three  curve  forms . In  the  some  way  the  square  deviation  S,! 
for  the  fitting  of  the  remaining  arrival  time  measurements  is  calculated 
for  the  three  expressions.  The  curve  form  g5.vi.ng  the  smallest  value  of 
S/S9*S"  is  the  best  one  by  this  criterion. 

This  has  been  done  for  each  shot.  For  equation  (5*2)  this  has  been  per- 
formed also  with  the  parameter  excluded.  (In  this  case  the  expression 

has  the  property  that  L = L when  t = t ).  The  results  when  D was  ex- 
eluded  and  when  was  included  turned  out  to  be  equally  good.  Each  gave 
the  smallest  value  of  S/S ,+S1'  in  about  50#  of  the  shots. 

In  30#  of  the  shots  equation  (5.2)  is  a better  fit  than  equation  (5.1) , 
and  in  90%  of  the  shots  equation  (5*2)  is  better  than  equation  (5*3). 
Thus,  equation  (5*2)  turns  out  to  fit  best  to  the  observed  arrival  times. 
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We  shall  in  the  following  sections  present  the  results  of  the  fitting  of 
equation  (5-2)  to  the  data  when  the  parameter  D is  included.  In  sections 
5-5  to  5.7  the  results  from  the  fitting  of  equation  (5.1)  to  the  data 
will  he  presented. 


5 . 2 Fitting  of  equation  (5,2)  to  the  measured  arrival  times 


A . B and  D have  teen  calculated  foi  each  shot,  and  listed  in  Tables  5.1a 
o ' o o 

and  0,  The  values  of  D and  S/(S:+SM)  are  also  listed  in  this  table. 

From  equation  (5,2)  the  shock  velocity  in  given  by  the  equation 


U 


dL 

at 


a 

o 


A 


B +(t“t  ) 
o o 


ar.d  the  acceleration  is 


(5.10 


dU 

at 


(B0+(t-t0))2 


(5-5) 


Ir.  position  L « 2C,  where  t * t , the  shock  velocity  exceeds  the  velocity 

o J 2 

of  sound  with  an  amount  A /B  , and  its  retardation  is  R - A /B  » When 

o c 000 

ve  shall  study  the  dependence  of  the  shock  on  scale  * charge  group  and 
roughness,  we  shall  take  Aq/B^  as  a measure  of  the  shock  strength  in  posi- 
tion N ~ 20,  and  DR^  to  represent  the  attenuation  of  the  shock. 


The  values  of  A /B  and  DR  are  calculated  for  each  shot.  Vie  shall  in 
00  o 

the  following  sections  take  the  values  of  A /B  and  DR  from  each  shot 

000 

ar.d  fit  simple  curve  forms  to  them,  in  order  to  study  their  qualitative 
Behaviour . 


5.3  Discussion  of  A /B 
’ — — — 


5.3,1  Dependence  of  A /B  on  q 

^00 

The  equation 


A 

_c 

B 


o 


to 
a e * 


(5.6) 


b and 


lit 


is  fitted  to  the  values  cf  A /B  > arid  the  resulting  values  of  a, 
the  correlation  coefficient  R are  listed  in  Table  5.2a  for  CFG.l  and 
Table  5.2b  for  CFG. 2.  b is  positive.,  showing  that  the  shock  strength 
in  position  II  = 20  increases  with  increasing  charge  weight,  b is  seen 
to  decrease  with  increasing  roughness.  This  is  consistent  with  what  was 
observed  in  section  H.2. 


Figure  5*1  shows  an  example  on  the  dependence  of  A /Bq  on  q. 

5.3.2  Dependence  of  A /B  on  rT 

o o L 


The  equation 


o 

B 

o 


a + brT 

JLi 


(5.7) 


is  fitted  to  the  values  of  A^/Bq » and  the  results  are  given  in  Table  5.3a 
for  CFC.l  and  in  Table  5. 3b  for  CFG. 2.  b is  seen  to  be  negative;  that  is, 
the  shock  strength  at  II  = 20  decreases  with  increasing  roughness.  This 
decrease  is  larger  for  the  larger  charge  groups. 


These  results  are  similar  to  those  obtained  in  section  U.2.  Gome  examples 
on  the  dependence  of  A /B  on  r are  shown  in  Figure  5*2. 

OO  J-4 


5.3.3  Dependence  of  A /B  on  D 

o o 


The  equation 


a + bD 


(5.8) 


is  fitted  to  the  values  of  A /B  . The  results  are  given  in  Table  5.^& 

0 0 

for  CFG.l  and  Table  5«^b  in  CFG, 2,  As  in  section  ’(.2,  a tendency  of  the 
shock  strength  at  K = 20  to  decrease  with  increasing  scale  is  observed, 
(t  < 0 in  29  of  33  cases).  This  tendency  seems  to  be  stronger  for  the 
larger  charge  groups.  It  is  difficult,  however,  in  this  case  to  see  any 
dependence  of  b on  r. . 


Figure  5 <3  shows  some  examples  on  the  deoendenc o A /B  on  D. 

o o 
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5. it  Discussion  of  DP. 

OT.mriT.  w • I.l  u ->  >M  -™~> . — Q 

, 2 . 

DR  53  D A /B  is  taken  to  represent  the  decay  in  the  shock  strength 
o o o 

as  the  Blast  propagates  ever  the  platform-  The  dependence  of  DR^  on 
charge  group ^ vail  roughness  in  the  tube  and  scale  will  nov  be  discussed 

5 *U.l  Dependence  of  DRq  on  q 

The  equation 

DR  ~ a eb<1  (5.9) 

o 

is  fitted  to  the  values  of  DR^*  The  results  are  listed  in  Table  5 * 5a  for 
CFG.l  and  Table  5. 5b  for  CFG. 2. 

In  section  4,3  the  attenuation  of  the  front  pressure  was  found  to  be 
nearly  independent  of  q.  From  Table  5*5  it  is  seen,  on  the  contrary, 

that  b > 0 in  14  of  18  cases,  indicating  that  the  decay  in  the  shock 

strength  increases  with  increasing  charge  group.  Hovcver,  no  systematic 
dependence  of  this  increase  on  vail  roughness  and  scale  is  observed. 

Figure  5.4  shows  an  example  on  the  dependence  of  DR^  on  q, 

5*4.2  Dependence  of  DR^  on 

The  results  from  fitting  the  equation 

DR  * a + brT  ( 5 .10) 

O b 

to  the  values  of  DR^  are  given  in  Table  5.6a  for  CFG.l  and  Table  5 -6b 
for  CFG. 2.  In  30  of  32  ea,ses  b is  found  to  be  negative  indicating  that 
the  decay  in  shock  strength  is  smaller  for  the  larger  values  of  vail 
roughness-  This  effect  was  also  observed  in  section  4-3. 

b is  seen  to  decrease  with  increasing  charge  group*  It  then  follows  that 
the  effect  of  smaller  attenuation  of  the  shock  with  increasing  roughness 
io  larger  for  larger  charge  groups* 

Figure  5*5  shows  some  examples  on  the  dependence  of  DRq  on  • 
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5.*+ i 3 Dependence  of  DR  on  D 

o 

The  equation 

DR  = a + bD  (5.11) 

o 

is  fitted  to  the  values  of  DF  . The  results  are  given  in  Table  5*7a  for 
CFG.l  and  Table  5,7b  for  CFG. 2.  In  30  of  33  cases  b is  negative.  This 
shows  that  the  decay  of  the  shock  strength  decreases  with  increasing 
scale.  The  similar  result  was  obtained  in  section  4.3. 

b seems  to  decrease  with  increasing  charge  group.  It  means  the  decrease 
of  the  attenuation  with  increasing  scale  is  larger  for  larger  charge 
groups.  The  tendency  here,  however,  is  weak. 

Some  examples  of  the  dependence  of  DRq  on  D are  shown  in  Figure  5*6. 


5 . 5 Fitting  of  equation  (5.1)  to  the  measured  arrival  times 

We  shall  also  briefly  present  the  results  from  the  fitting  of  equation 
(5.1)  to  the  measured  arrival  times,  t in  equation  (5.1)  is  taken  to 
be  the  measured  arrival  time  in  position  N --  20.  From  equation  (5.1)  the 
velocity  of  the  shock  front  is  given  by 

U = D(V2$2(K~20>  + 363(H-20)2r1  (5.12) 

and  its  acceleration  is 


au 

dt 


-D(B,+2B2(H-20) 


333(N-20)2)~3(2g2+6&  (JJ-20))  (5.13) 


The  shock  velocity  at  N « 20  is  then 


U 

o 


D 

e7 


and  the  retardation  is 


R 


o 


+ 


3 


(5.1*0 


(5.15) 


IT 


Since  our  main  interest  now  is  to  compare  the  results  from  the  fitting 

of  equation  (5.1)  to  the  measured  arrival  times  with  the  results  from 

the  fitting  of  equation  (5*2),  we  shall  concentrate  the  discussion  on 

U and  the  "scaled"  retardation  in  N = 20.  DR  . It  must  be  borne  in 
o ' o 

mind,  however,  that  there  is  an  additional  parameter  3 . 


The  values  of  U , DR  and  3~  are  given  in  Table  5 >8?-  for  CFG.l  and. 
c>  o 3 

Table  5.8b  for  CFG, 2. 


5.6  Discussion  of  U 

Q 

5 .6.1  Dependence  of  U on  q 
The  equation 

Uo  = a e’Dq  (5.16) 

is  fitted  to  the  values  of  U^,  and  the  results  are  listed  in  Table  5»9a 
for  CFG.l  and  Table  5. 9b  for  CFG. 2.  As  in  section  5*3.1  b is  found  to 
be  positive,  and  it  has  a tendency  to  decrease  with  increasing  wall  rough 
ness  and  scale. 

5.6.2  Dependence  of  U on  rT 

o L 

The  equation 

Uq  « a + brL  (5.17) 

is  fitted  to  the  values  of  U and  the  results  arc  given  in  Tables  5 -1C. 
b is  found  to  be  negative,  and  it  has  larger  negative  values  for  larger 
charge  groups . 

The  results  are  similar  to  those  observed  in  section  5*3.2. 

5.6.3  Dependence  of  lTo  on  D 

The  results  from  fitting  the  equation 


a + ID 


(5.18) 
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to  the  values  of  arc  listed  in  Tables  5 .11.  b is  negative  in  23  of 
33  cases,  so  the  tendency  of  the  shock  strength  in  N ~ 20  to  decrease 
with  increasing  scale  is  observed  also  in  this  case.  As  in  section  5.3.3 
this  tendency  seems  to  be  stronger  for  larger  charge  groups. 


5.7  Discussion  of  DR 
' o 

5.7.1  Dependence  of  DRq  on  q 

Equation  (5.S; ) is  fitted,  to  the  values  of  DRq  and  the  results  are  listed 
in  Tables  5.12.  As  in  section  5*^.1  a weak  tendency  of  DRq  to  increase 
with  increasing  charge  group  is  observed.  It  is  difficult  to  see  any 
systematic  dependence  of  this  increase  on  scale  and  roughness. 

5.7.2  Dependence  of  DR  on  r 

O L* 

Equation  (5*10)  is  fitted  to  the;  values  of  DR^,  and  uhe  results  are 
listed  in  Tables  5-13.  DR  is  found  to  decrease  with  increasing  rT , just 

O Jj 

as  in  section  5.^.2.  A weak  tendency  of  this  decrease  to  become  larger 
for  larger  charge  weights  is  also  observed,  while  it  is  difficult  to  find 
any  dependence  on  D. 

5.7*3  Dependence  of  DF^  on  D 

At  last,  equation  (5 .11)  is  fitted  to  the  values  of  DRq,  and  from  the 
results,  listed  in  Tables  5.1^,  DR  does  not,  as  in  section  5*^.3,  show 
any  tendency  to  decrease  with  increasing  scale.  This  difference  is  due 
to  the  additional  fitting  parameter  in  Equation  (5.1).  When  the  scale- 
dependence  of  this  parameter  is  taken  into  account , the  retardation  of 
the  shock  front , and  hence  the  attenuation  of  the  front  pressure,  is 
found  to  decrease  with  increasing  scale. 

5 . 6 Summary 

In  this  and  the  proceeding  chapter  several  curve  fitting  procedures  were 
performed  to  interpret  the  arrival  time  measurements , They  give  all  the 
same  qualitative  results.  This  assures  that  the  conclusions  are  not  only 
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resulting  from  the  curve  fitting  procedures  used,  but  reveal  the  physical 
content  of  the  data. 

The  main  results  can  be  summarized  as  fellows:  At  the  beginning  of  the 
platform  (N  = 20),  the  front  pressure  increases  with  increasing  charge 
group,  and  decreases  with  increasing  scale  and  roughness.  The  attenuation 
of  the  front  pressure,  however,  also  decreases  with  increasing  scale  and 
roughness . 


6 


THE  90  mbar  LIMIT 


6 . 1 Determination  of  the  90  mbar  limit 

It  is  of  interest  to  find  the  distance  from  the  tube  outlet  in  which  the 
front  pressure  falls  below  a certain  value,  for  example  50  mbar.  Since 
equation  (5.2)  gives  the  best  fit  to  the  measured  arrival  times,  this 
will  be  used  to  determine  the  50  mbar  limit. 

By  use  of  equation  (5*2)  and  the  Rankine-Hugenoit  equation,  the  50  mbar 
limit  is  found  to  be  expressed  as 


x = 20  + ~~{Vr.6  Ao-aoBc+Ao  In (1*7 


A 


(6,1) 


x measures  the  distance  in  number  of  tube  diameters.  The  values  of  x 
obtained  by  Equation  (6.1)  are  listed  in  Table  6.1  for  the  different 
charge  groups,  degrees  of  wall  roughness  and  scales. 


6.2  Discussion  of  x 


6.2,1  Dependence  of  x on  q 


The  equation 


b 


x = a q 


is  fitted  to  the  values  of  x,  and  the  results  arc  given  in  Table  6.2a  for 
CFG.l  and  Table  6,2b  for  CFG, 2.  All  the  values  of  b are  found  to  be  be- 
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tveen  0 and  3.,  i.  e,  about  O.k  for  CFG . 1 and  0.5  for  CFG.2S  shoving  that 
x increases  with  increasing  q.  Figure  6.1  shows  an  example  on  the 
dependence  of  x on  q.. 


It  is  difficult  to  find  any  dependence  of  b on  D and  rL< 


L 


6.2.2  Dependence  of  x on  r7 
The  equation 

x = a + hr. 


(6.3) 


is  fitted  to  the  values  of  x,  and  the  results  are  given  in  Table  6.3. 
b is  found  to  be  negative  in  26  of  33  cases,  indicating  that  x decreases 
with  increasing  roughness  in  the  tube.  The  decrease  seems  to  be  larger 
for  the  larger  charge  groups. 

As  discussed  in  the  proceeding  chapters,  the  effect  of  the  wall  rough- 
ness in  the  tube  is  to  flatten  the  pressure-time  profile  of  the  blast. 

This  leads  to  smaller  values  of  the  front  pressure  at  the  beginning  of 
the  platform  in  the  cases  of  large  roughness,  but  also  to  smaller  attenua- 
tion of  the  front  pressure  as  the  shock  propagates  over  the  platform. 

For  the  50  mbar  limit  the  effect  of  the  roughness  still  is  to  give  smal- 
ler values  of  x with  larger  degrees  of  rovighness. 

In  Figure  6.2  are  shown  some  examples  on  the  dependence  of  x on  r ^ . 


6.2.3  Dependence  of  x on  D 


The  equation 


x 


a + bD 


(6J0 


is  fitted  to  the  values  of  x.  The  results  are  listed  in  Tables  60. 


b is  positive  in  26  of  33  cases,  showing  that  the  distance  to  the  point 
where  the  pressure  falls  below  50  mbar  increases  with  increasing  scale, 
b > 0 in  2k  of  33  cases  is  a relative  weak  tendency.  This  is  probably 
due  to  the  effect  of  the  tube  outlet.  The  front  pressure  in  position 
N = 20  on  the  platform  was  found  to  decrease  with  increasing  scale.  The 
attenuation  of  the  front  pressure  as  the  shock  propagates  over  the  plat- 
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form,  was,  however,  also  found  to  decrease  with  increasing  scale,  so  than 
for  larger  values  of  N increasing  front  pressure  with  increasing  scale 
is  observed,.  The  value  cf  N for  which  this  occurs,  will  depend  on  the 
pres sure “time  profile  of  the  blast,  but  no  systematic  analysis  of  this 
has  been  performed. 

Figure  6,3  shows  some  examples  on  the  dependence  of  x on  D. 


7 COMPARISON  WITH  FULL  SCALE  TESTS 


We  shall  briefly  compare  the  results  obtained  in  Norwegian  full  scale 
tests  at  Raufoss  and  Storfosen  with  the  present  model  tests.  The  Raufoss 
and  Storfosen  experiments  are  described  in  (6)  and  (7)?  respectively. 

There  are  large  differences  both  in  charge  groups  used  and  in  geometry 
between  these  full  scale  tests  and  the  model  experiments.  This  makes  a 
direct  comparison  difficult,  Except  for  one  shot , where  a direct  compari- 
son is  performed,  ve  have  therefore  only  taken  the  front  pressures  ob- 
tained in  the  Raufoss-  and  Storfosen  shots  and  fitted  equation  (^.l)  to 
them.  The  results  are  given  in  Table  7*1*  The  curve  fitting  is  based  on 
measurements  in  2,  ^ or  5 positions  outside  the  tunnel  outlet,  N being 
between  10  and  80« 

The  tunnel  diameter  at  Raufoss  is  D ~ 2.65  rti,  corresponding  to  scale 
53:100.  At  Storfosen  2 tunnels  were  used.  Both  tunnels  had  bends,  and  in 
both  tunnels  the  diameter  changed.  Ve  shall  associate  a diameter  D = 

2.85  m to  the  one  and  a diameter  D » 1*  m to  the  other,  corresponding  to 
scales  57:100  and  CC:100,  respectively. 

In  the  diacussioh  in  section  U.3  the  attenuation  B was  found  to  decrease 
with  increasing  scale.  By  comparing  the  values  of  B from  the  full  scale 
case  in  Table  7*1  with  the  values  of  B from  the  model  tests  in  Tables  H.i, 
the  attenuation  in  the  full  scale  tests  is  of  the  same  order  of  magni- 
tude or  somewhat  smaller.  This  may  be  due  to  the  different  values  of 
charge  group  q,  although  in  section  U.3.  B was  found  to  be  nearly  in- 
dependent of  q. 

from  Table  7 I with  the  values  of  P from 

c 


By  comparing  the  values  of  P 

o 
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Table  4.1,  P is  found  to  be  smaller  in  the  full  scale  tests.  This  is 
o 

reasonable,  since  the  charge  groups  are  smaller  in  the  full  scale  tests. 

As  mentioned  there  are  several  important  differences  between  these  full 

scale  experiments  and  the  model  tests,  and  this  makes  a direct  comparison 

difficult.  It  is  shot  no  9 from  Raufoss  which  is  best  suited  for  a direct 

comparison,  and  even  in  this  case  the  differences  are  large.  The  charge 

weight  in  this  shot  is  Q » 1C00  kg,  corresponding  to  charge  group 

q ss  53.7  kg/rf"* . The  pressure  as  a function  of  distance  from  the  tunnel 

outlet  is  shown  in  Figure  7-l>  The  Raufoss  tunnel  is  closed  in  one  end. 

as  our  CFG. 2,  and  the  charge  is  detonated  in  the  tunnel,  as  in  our  CFG.l. 

The  Raufoss  shot  should  then  be  compared  with  shots  with  charge  group 

2x53.7  kg/m  ~ .107  kg/m  in  CFG.l  in  our  model  tests.  We  have  not  used 

exactly  this  charge  group,  but  a comparison  can  be  performed  by  inter- 

3 3 

polation  from  the  shots  with  q - 80  kg/m  and  q =:  loO  kg/m  . The  result- 
ing pressures  as  function  of  distance  from  tube  outlet  for  scales  1:100, 
2.100  and  3:100  are  also  shown  in  Figure  7*1*  It  is  assumed  that  rT  ~ 

0.06  is  the  wall  roughness  which  best  corresponds  to  the  roughness  in 
the  Raufoss  tunnel. 

As  is  seen  in  Figure  7.1,  that  while  the  pressure  increases  with  increas- 
ing scale  in  the  model  tests,  the  pressure  is  much  lower  in  the  Raufoss 
shot . 

It  seems  reasonable  to  explain  this  deviation  between  the  model  tests 
and  the  Raufoss  shot  by  the  larger  differences  between  the  full  scale 
case  and  the  models,  of  which  the  most  important  are: 

“ A possibly  larger  large-scale  and  small-scale  roughness  in  the  rock 
tunnel  at  Raufoss 

- Vegetation  and  irregularities  in  the  terrain  might  be  of  greater 
importance  than  irregularities  on  the  platform 

~ There  is  a side-chamber  and  a bend  in  the  Raufoss  tunnel 

- Energy  in  the  blast  wave  is  used  fox'  evaporation  of  water  in  the  tunnel 

- In  the  models  there  are  a vertical  steel  plate  above  the  tube  outlet,  . 
while  at  Raufoss  the  hill  above  the  tunnel  is  not  very  steep 

- The  tunnel  length  at  Raufoss  is  only  3^  tunnel  diameters;  the  effect 
of  this,  however,  should,  be  in  the  opposite  direction 
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8 CONCLUSIONS 


For  the  measurements  in  the  tube  section  a deviation  from  the  scaling 
lavs  vas  observed  (l).  It  was  seen  that  the  front  pressure  increased 
with  increasing  scale,  mainly  due  to  the  fact  that  the  attenuation  of 
the  front  pressure  decreased  with  increasing  scale. 

In  chapter  3 of  the  present  report  the  tube  outlet  was  seen  to  have  the 
effect  of  reversing  this  deviation;  at  the  first  measurement  stations 
after  the  tube  outlet  the  fro?it  pressure  decreased  with  increasing  scale. 
In  other  words,  there  is  a larger  pressure  drop  in  the  larger  scales, 
when  the  shock  propagation  changes  from  one-dimensional  in  the  tube  to 
three-dimensional  on  the  platform.  This  effect  was  seen  to  depend  on  the 
wall  roughness  in  the  tube,  indicating  that  the  pressure-time  profile 
of  the  blast  is  of  importance  for  the  magnitude  of  the  effect. 

The  effect  does  not  seen  to  result  from  a scale  dependent  loss  of  energy 
in  the  tube  outlet,  but  rather  from  a scale  dependent  redistribution  of 
energy.  The  attenuation  of  the  front  pressure  as  the  shock  propagates 
over  the  platform  is  smaller  in  the  larger  scales,  and  at  sufficiently 
large  distances  from  the  tube  outlet  the  front  pressure  increases  with 
increasing  scale.  This  effect  makes  it  difficult  to  extrapolate  the 
scale  dependence  of  the  50  mbar-limit,  for  example,  to  large  scales. 

In  chapters  k and  5 where  the  dependence  of  the  attenuation  of  the  front 
pressure  and  the  front  pressure  in  position  N = 20  on  scale,  charge  group 
and  roughness  was  studied,  the  attenuation  vas  found  to  depend  on  the 
roughness  in  the  tube.  This  again  shows  the  importance  of  the  whole  pres- 
sure-time history  for  describing  the  blast. 

The  time  frame  of  this  project  has  not  made  a study  of  the  pressure-time 
profiles  possible,  although  they  are  recorded* 

In  summary , a relative  complicated  dependence  of  the  shock  strength  cn 
scale  is  observed.  This  supports  the  conclusion  in  (l),  that  the  simple 
scaling  lav/s  must  be  used  with  great  care . 
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Figure  3* la  Pressure  drop  in  the  tube  outlet 


Pressure  drop  -in  the  tube  outlet 


Figure  3.1b 
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Figure  5.1 


A-./-V v as  function  of  charge-  rrouu , 

('xhe  curve  is  the  result  "from  fitting  eq  (5,6)  to  the 

values  of  A /E  ) 
c'  o 
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Figure  5.2  A/r/B^  as  function  of  vail  roughness,  rj 

(The  curves  are  the  result  from  fitting  eq.  (5.7)  to  the 

values  of  A /B  ) 
o o' 


Figure  5.2 


A- /Bc  as  fuac tion  of  tube  diameter,  D 

(The  curves  are  the  result  from  fitting  eq  (5.8)  to  the 

values  of  Aq/Bq) 
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Fi ^ur  e 5 - 5 RnD  as  a function  of  wall  roughness,  i y( 

(The  curves  are  the  result  from  fitting  eq  (5.10)  to 
the  values  of  RqD) 
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Figure  5.6  R0D__ as  functi  or.  of  tubs  diameter  D 

(The  curves  are  the  result  from  fitting  eq  (5*11)  to  the 
values  of  RqD) 
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Figure  6.2  Dependence  of  the  50  nbar  lirnit  on  vail  roughness,  rj 

(The  curves  are  the  result  from  fitting  eq.  (6.3)  to  the 
values  of  x) 
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Figur e 6.3  Dependence  of  the  50  rr/har  Unit  on  tute  aiaraet_er^  JD 

(The  curves  arc  the  result  fremt  fitting  eq  (6*H)  to  the 
values  of  x) 
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Table  5.8a  Results  from  fitting  of  eg  (5.1)  to  the 
measured  arrival  times 

0vRo  is  multiplied  by  1G~*  and  go  bv  lo£ 


51* 


- 

CFG.2 

0 - O.Ot 

D = 0.18 
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Table  6.1  The  fifty  millibar  limit 


62 


CFG.?  X*aqb 
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No  of  measure- 
ment positions 
included  in 
curve  fitting 

P = Po 

W 

D(m) 

SHALE 

Q 

G 

«is  — 

Po 

B 

R 

STORFGSEN 

4.00 

80:100 

72.4 

1.13 

4 

0.034 

1.355 

0.94 

STORFOSEN 

4.00 

80:130 

375.3 

0.86 

5 

0.C97 

1.434 

0.38 

STORFOSEN 

2.85 

57:100 

151.3 

G.&9 

2 

C.033 

1.293 

- 

STORFOSEN 

2.85 

57:103 

382.5 

127)7 

*\ 

C 

0,053 

1.128 

- 

STORFOSEN 

2.85 

57:100 

1490 

64.1 

2 

0.140 

G.Q93 

- 

RAUF OSS 

2.65 

53:160 

130 

5.37 

4 

C.C33 

6.984 

1.00 

RAUFOSS 

2.65 

52:1  GH 

300 

16.12 

1 

0.046 

1.073 

- 

R4UF0SS  * 

2.65 

53:100 

1 0lfO 

53.7 

3 

0.090 

t.315 

0.95 

R4UF0SS  * 

2.65 

53.100 

1000 

53.7 

3 

0.065 

1.120 

0.9S 

* The  charges  are  placed  at  different  distances  from  the  closed  end  of  the  tunnel. 


Table  7*1  Results  from  sor^e  Nerve ri an  fall  scale  tests.. 


i 


f 


